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Abstract 

The purpose of this paper is to show that Non- Archimedean Mathe- 
matics (NAM), namely mathematics which uses infinite and infinitesimal 
numbers, is useful to model some Physical problems which cannot be de- 
scribed by the usual mathematics. The problem which we will consider 
here is the minimization of the functional 

E(u,q) = i JjVu(x)\ 2 dx + u(q). 

If n C ~R N is a bounded open set and u G Co(f2), this problem has no 
solution since inf E(u, q) = — oo. On the contrary, as we will show, this 
problem is well posed in a suitable non- Archimedean frame. More pre- 
cisely, we apply the general ideas of NAM and some of the techniques of 
Non Standard Analysis to a new notion of generalized functions, called 
ultrafunctions, which are a particular class of functions based on a Non- 
Archimedean field. In this class of functions, the above problem is well 
posed and it has a solution. 

Mathematics subject classification: 26E30, 26E35, 35D99, 35J57. 

Keywords. Non Archimedean Mathematics, Non Standard Analysis, 
ultrafunctions, Delta function, Dirichlet problem. 
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1 Introduction 

By Non-Archimedean Mathematics (NAM) we mean mathematics based on 
non-Archimedean fields, namely ordered fields which contain infinite and in- 
finitesimal numbers. We believe that, in many circumstances, NAM allows to 
construct models of the physical world in a very elegant and simple way. In the 
years around 1900, NAM was investigated by prominent mathematicians such 
as Du Bois-Reymond [6], Veronese [12], David Hilbert j9] and Tullio Levi-Civita 
[5], but then it has been forgotten until the '60s when Abraham Robinson pre- 
sented his Non Standard Analysis (NSA) [TT]. We refer to Ehrlich [7] for a 
historical analysis of these facts and to Keisler [TU] for a very clear exposition 
of NSA (see also Q], @]) 

The purpose of this paper is to show that NAM is useful to model some 
Physical problems which cannot be described by the usual mathematics even if 
they are relatively simple. 

The notion of material point is a basic tool in Mathematical Physics since 
the times of Euler who introduced it. Even if material points do not exist, 
nevertheless they are very useful in the description of nature and they semplify 
the models so that they can be treated by mathematical tools. However, as 
new notions entered in Physics (such as the notion of field), the use of material 
points led to situations which required new mathematics. For example, in order 
to describe the electric field generated by a charged point, we need the notion 
of Dirac measure 5 q , namely this field satifies the following equation: 



where A is the Laplace operator. 

In this paper, we will describe a simple problem whose modelization requires 
NAM. Let f2 C R 2 be an open bounded set which represents a (ideal) membrane. 
Suppose that in f2 is placed a material point P, which is left free to move. 

Suppose that the point has a unit weight and the only forces acting on it 
are the gravitational force and the reaction of the membrane. If q 6 f2 is the 
position of the point and u(x) represents the profile of the membrane, it follows 
that equation ([1]) holds in f2 with boundary condition u — on <9f2. 

The question is: which is the point qo <E Q that the particle will occupy? 
The natural way to approach this problem would be the following: for every q e 
fi, the energy of the system is given by the elastic energy plus the gravitational 
energy, namely 



Alt = S q 



(1) 
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If the couple (uo,go) minimizes E, then go is the equilibrium point. For every 
ggil, let u q (x) be the configuration when P is placed in q 1 namely the solution 
of equation ([!]). So the equilibrium point qo is the the point in which the function 



has a minimum. 

In the classical context, this "natural" approach cannot be applied; in fact 
u q (x) has a singularity at the point q which makes u(q) not well defined and 
the integral in ^ to diverge. On the contrary, this problem can be treated in 
NAM as we will show. In fact, since infinite numbers are allowed, we will be 
able to find a minimum configuration for the energy 

In order to pursue this program, we apply the general ideas of NAM and 
some of the techniques of NSA to a new notion of generalized functions which 
we have called ultrafunctions (see 0). Ultrafunctions are a particular class 
of functions based on a superreal field R* D R. More exactly, to any continuous 
function / : M. N — > R we associate in a canonical way an ultrafunction /$ : 
(R*) — > R* which extends /; but the ultrafunctions are much more than the 
functions and among them we can find solutions of functional equations such as 
equation |T]) which are defined in every point of Q* C (R*)^ . Thus, the energy 
(|2j) is well defined for every ultrafunction even if it might assume infinite values. 

Now we itemize some of the peculiar properties of the ultrafunctions: 

• the space of ultrafunctions is larger than the space of distributions, namely, 
to every distribution T, we can associate in a canonical way an ultrafunc- 
tion T$ (for details see [2]); in particular the Dirac measure can be repre- 
sented by an ultrafunction S q (x) and, for every ultrafunction u, we have 
that 



• similarly to the distributions, the ultrafunctions are motivated by the need 
of having generalized solutions and also by the need to model extreme 
physical situations which cannot be described by functions defined in 18L N ; 
however, while the distributions are no longer functions, the ultrafunctions 
are still functions even if they have larger domain and range; 

• unlike the distributions, the space of ultrafunctions is suitable for non 
linear problems such as the one described above; 

• if a problem has a unique classical solution u, then u$ is the only solution 
in the space of ultrafunctions; 

• the main strategy to prove the existence of generalized solutions in the 
space of ultrafunction is relatively simple; it is just a variant of the Faedo- 
Galerkin method. 

Before concluding the introduction, we refer to [3] and to [5] where other 
situations which require NAM are presented. 



F(q) = E(u q ,q) 



(3) 
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1.1 Notations 

Let fl be a subset of M. N : then 

• C (£1) denotes the set of real continuous functions defined on ft; 

• Cq (jTj denotes the set of real continuous functions on fl which vanish on 

on-, 

• C k (f2) denotes the set of functions defined on ft c which have contin- 
uous derivatives up to the order fc; 

• C$ (0) = C k (IT) n C (H) ; 

• T> (O) denotes the set of the infinitely differentiable functions with compact 
support defined on 17 C M. N ; 

• H 1 ^!) is the usual Sobolev space defined as the set of functions aeL 2 (f2) 
such that Vu G L 2 (fi) ; 

• H^(Q) is the closure of V (Q) in H^Q); 

• iJ — 1 (O) is the topological dual of if^fi). 

2 The ultrafunctions 

In this section we briefly recall the notion of A-limit and of ultrafunction which 
have been introduced in [2J. 

2.1 The A-limit 

The idea behind the concept of A-limit is the following: let U denote a "math- 
ematical universe" (which will be precisely introduced in definition ((3])), and T 
the set of finite subsets of U, ordered by inclusion. The A-limit can be thought 
as a way to associate to every net ip : J- — > K a limit lim^A) that satisfies a 

few properties of coherence. 

These limits will be elements of a Non- Archimedean field K; since this leads to 
work in such fields, we recall a few basic facts and definitions: 

Definition 1 Let K be an ordered field. Let (6t We say that: 

• £ is infinitesimal if for all n S N |£| < \ ', 

• £ is finite if there exists neN such as |£| < n; 

• £ is infinite if for all n £N |£| > n. 

Definition 2 An ordered field K is called non- Archimedean if it contains an 
infinitesimal £ 0. 
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We are interested in fields that extend R: 



Definition 3 A superreal field is an ordered field K that properly extends R. 

Since R is complete, it is easily seen that every superreal field contains in- 
finitesimal and infinite numbers. 

In order to precise the notion of A- limit, we need to define the notion of "math- 
ematical universe" . For our applications, we take as mathematical universe the 
superstructure on R: 

Definition 4 The superstructure on R is 

oo 

U = |J U„ 

71=0 

where U n is defined by induction as follows: 

U = R; 
U„+i = V n UP(V n ) 

Here V (E) denotes the power set of E. If we identify the couples with 
the Kuratowski pairs and the functions and the relations with their graphs, U 
formalizes the intuitive idea of mathematical universe. 

We denote by T the set of finite subsets of U. Ordered with the relation of 
inclusion, T becomes a direct set; following the usual nomenclature, we call net 
(with values in E) any function ip : T — > E. 
Following [5] , we introduce the A- limit axiomatically: 

• (A-l) Existence Axiom. There is a superreal field K Z> R such that for 
every net ip : J- — > R there exists a unique element L £ K called the 
"A-limit" of ip. The A-limit will be denoted by 

L = lim (p(\). 

Moreover we assume that every £ € K is the A-limit of some net ip : J- — > 



• (A-2) Real numbers axiom. // ip(X) is eventually constant, namely 
3Ao £ J: VA D Ao, y(A) = r, then 

lim tp(\) = r. 

AfU 

• (A-3) Sum and product Axiom. For all <p, %p : T — > M: 

Atu^ + AtS^ = 1™(^ A ) + V>(A)); 
lWA)-lim^(A) = lim(^(A)^(A)). 

AfU AtU AfU 
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Theorem 5 The axioms (K-l),(K-2),(K-3) are consistent. 

Proof. This is the content of Theorem 7 in [5] . 

□ 

We say that a net ip : J- — > U is bounded if 

3n G N such that, VA G J 7 , <p(\) G U„. 

The notion of A-limit can be extended to bounded nets by induction on n: 
for n = 0, lim 93(A) is defined by the axioms (A-l),(A-2),(A-3); so by induction 

we may assume that the limit is defined for n — 1 and we define it for a net 
(p : T — > U„ as follows: 

limw(A) = llim^(A) I ip : J" ^ U„_i and, VA G A, t£(A) G w(A) 

A set that is a A-limit of sets is called internal. The A-limit provides a way to 
extend subset of M and functions defined on (subsets of) K to K: 

Definition 6 Given a set E C K let ce '■ T — >• U &e i/ie nei smc/i </ia< VA 6 J" 
ce(A) = £. TTierc 

:= limc B (A) - (w(A) | ^(A) G E 

AfU A^U 

is called natural extension of E. 

Using the above definition we have that 

K = R*. 

A function / can be extended by identifying / and its graph, and this extension 
satisfies the following properties: 

Theorem 7 For every sets A,BeV, the natural extension of a function 

f-A^B 

is a function 

f* : A* ->• B*\ 
moreover for every ip : A n V (A) — > A, we have that 



lirn f(<p{\)) = f* ( lim^(A) 

AfU V^AtU 

A property that is natural to ask for the A-limit of a net tp is that some 
properties of the limit can be deduced from the properties of p. This is ensured 
by the following important theorem: 



G 



Theorem 8 (Leibnitz Principle) Let 1Z be a relation in U„ for some n > 
let <p,ip : T ->■ U n . // 



VA g J 7 , <^(A)^(A) 

The last key concept that we need is that of hyperfinite set: 

Definition 9 An internal set is called hyperfinite if it is the A-limit of finite 
sets. 

All the internal finite sets are hyperfinite, but there are hyperfinite sets which 
are not finite, e.g. the set 

M° := lim(R n A) 

AfU 

is not finite. The hyperfinite sets are very important since, by Leibnitz Prin- 
ciple, they inherit many properties of finite sets; e.g., R° has a maximum and 
a minimum clement, and every internal function (i.e. a function such that its 
graph is an internal set) 

/ : R° -> K* 

has a maximum and a minimum as well. Intuitively, hyperfinite sets can be 
thought as having an hyperfinite number /? of elements, where (3 is an element 
ofN*. 

Given a set A G U we denote by A° its hyperfinite extension: 

A° = lim(A n A). 

atu 

By this construction, if a hyperfinite set consists of numbers, or vectors, it 
is possible to add all its elements. Let 

A := YnnA\ 

be a hyperfinite set; the hyperfinite sum of the elements of A is defined as 
follows: 

Ea = lim >^ a 
AtU ^ 

aeA agAj, 

In particular, if A = {a\, ap} consists of (3 elements, with /? G N*, we use the 
notation 

a 

aeA j=l 
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2.2 Definition of the ultrafunctions 

Let f2 be a subset of R N , and let Vq (ft) be a vector space such that V(Q) C 
v G (n) CC(H)nL 2 (0). 

Let </3y G (o) be the net such that, for every Ae J, <Pv a (ii)W = Vx(^)> where 

V X {Q) = Span(V G {Q)r)X). 
Definition 10 The set of ultrafunctions generated by Vc(f2) is 
V(fi) = timVx(n) = 5pon(Vfe(n)°); 

any element u(x) of V(Ct) is called ultrafunction and Vq is called the 
generating space. 

Observe that, being the A-limit of a net of vectorial spaces of finite dimen- 
sions, V(Q) is a vectorial space of hyperfinite dimension. Its dimension, that we 
denote by /?, is 

P = \imdim(V x (Cl)). 

The ultrafunctions are A- limits of continuous functions in V\(Q), so they are 
internal functions 

u ■. n* -> C*. 

(we recall that a function is called "internal" if it is a A-limit of functions). 

Notice that V(fl) inherits an Euclidean structure that is the A-limit of the 
Euclidean structure of every space V\(d) given by the usual L 2 (f2) scalar prod- 
uct; also, since V(Q) is a subset of L 2 (Q)* , it can be equipped with the following 
scalar product 



(u,v) — / u(x)v(x) dx. 
Jn 

where J Q is the natural extension of the Lebesgue integral considered as a func- 
tional. 

Being a vectorial space of hyperfinite dimension, V(Q.) admits an hyperfinite 
orthonormal basis {ei(x) \ i < P}. Having fixed a basis, we can make two im- 
portant constructions in an explicit form. The first is the extension to V(Sl) of 
continuous functions f(x) such that 

\/v(x) e V(Q) - oo < / f*(x)v(x)dx < +oo. 
Jn 

Let f(x) be such a function, and let $ denote the orthogonal projection 

$ : C(fi)* -> V(Cl). 
We call canonical extension of f(x) the ultrafunction 
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U{x) = <&(/*(*)). 

Observe that = /* <^> f(x) 6 Vg(S1) as expected, and that for every 
function f(x) the following holds: 

f*(x)v(x)dx = / /$(^)u(a;)da;. 

In terms of the basis {ei(x) \ i < /?}, the operator $ has the following 
expression: 

*(/(«)) - /*(*) = E ( f no^m) (4) 

The second important construction regards the Dirac delta functions: 

Theorem 11 Given a point q G S7, there exists a unique function S q (x) in V(f2) 
such that 



□ 



Vv e V(Q), j S q {x)v(x) dx = v(q). (5) 
Proof. The proof can be found in [2], Theorem 23. 



is called the Dirac ultrafunction in V(Sl) concentrated in q. In terms 
of the basis {ei(x) \ i < /3}, the 5 q has the following expression: 



S q( x ) = ^2 e t(v) e i( x ), (6) 
i=l 

which validity can be checked with a direct calculation. 

Remark 12 We observe that, in the context of ultraf unctions, the Dirac ultra- 
functions are actual functions, while in the classical theory of functions they are 
distributions. For example, in the ultrafunction context it makes perfect sense 
to consider objects like 5 q (x) 2 , 6 q (x) — 1, 6 q (x) ■ 8 q '(x) and so on. 

3 The model problem 

In this section we want to solve the problem described in the introduction via a 
"natural" approach that can not be applied in the classical framework, while it 
can be applied in the ultrafunction setting. We begin by describing the Dirichlet 
problems in the framework of ultrafunctions. 
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3.1 The Dirichlet problem 

Let fl be a bounded open set in M. N , and consider the Dirichlet problem: 

-An = f(x) for ieJ] ^ ' 

When dfl and f(x) are smooth problem ([7]) has a unique solution. Otherwise, 
in the classical Sobolev approach, problem ([7]) is trasformed in the following: 

u g ^(O) 

-Au = f(x) w 

Problem (|5J has a unique solution whenever f2 is a bounded open set and 
f(x) is in in this case the equation — Au = / is required to be satisfied 

in a weak sense: 



/ uAip dx — / ftp dx \/ip G T>(il) 
Jn Jo. 



Also, the solution u(x) given by this procedure is not a function but an 
equivalence class of functions defined a.e. in Q. 

In the approach with ultrafunctions let Vj 3 2 (0) be the space of ultrafunctions 
generated by Cq(Q). Problem ([7]) can be rewritten as follows: 

u G V Q 2 (il) 

-A 9 u = f(x) for x G n* { ' 

where A$ = $ o A* : V^(Cl) ->• F 2 (ft). 

Observe that now we are solving the problem in an hyperfinite space, and by 
Leibnitz Principle it follows that there is an unique solution for every f(x) G 
Vq (fi) (for the details see [2] , Theorem 27) . The idea of the proof is the following. 
The solution can be constructed by first finding a solution U\(x) in each finite 
dimensional space (Vq(£1))\ = Span(C 2 (Vl) H A), and then taking the A-limit 

u(x) = ]hau\(x). 

The solution u{x) is an ultrafunction defined for every x G Q* and, since 

Vx G 90, VA eJnC 2 (fi), u\{x) = 0, 

it follows by Leibniz principle that Vx G 9(0*), u(x) = 0. 

So satisfies the pointwise boundary condition, a result that is not true 
in the Sobolev approach. Finally, when problem ([Jj has a solution s(x) G C 2 (0), 
then 

— s* (x) 

and, when problem ([8]) has a solution g(x) G -ffg(O), then we have that 
g(x)v (x) dx ~ / u(x)w(a;) <±e Vw(.t) G Co(O) 
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3.2 A solution by mean of ultrafunction 

Now let us consider a minimization problem inspired by the one which we have 
discussed in the introduction. Let ft C R N be an open bounded set; we want 
to find a function u defined in il (with u = on 90) and a point geSJ which 
minimize the functional 



E(u,q) = i Jj\7u{x)\ 2 dx + u(q) 



It is well known that this problem has no solution in Cq (ft) and it makes no 
sense in the space of distributions. On the contrary it is well defined and it has 
a solution in Vg 2 (f2). More exactly, we have the following result: 

Theorem 13 For every point q G £1* , the Dirichlet problem 

A$m = 5 q for x G O* 
u(x) = for x G <90* 

/ias a unique solution u q G V" 2 (O) whose energy E(u q ,q) G R* is an infinite 
number; moreover there exists qo G 0* s«c/i i/iat 

E(u qo ,q Q )= mmE(u q ,q) = min E(u,q). 

gen* gen*, Mev 2 (n) 

Proof: First of all we observe that 

min £7(u g , q) = _ min E(u,q) 
gen* qesi*, uev^(n) 

since every stationary point (u,q) of E(u,q) satisfies — A$(u) = 

To minimize E(u, q) we use the Feado-Galerkin method, namely the finite 
dimensional reduction. First of all, for every A 6 J, we solve the following 
problem in (V 2 {Q)) X = Span(C 2 (n) n A): 

ue(V a 2 (Q)) x 
J Au v dx — J S q v dx for every v G (Iq 2 (0))a 

This problem has a unique solution u q .\(x) for every A G T fl C 2 (0), since 
(V r 2 (0))A is a nonempty finite-dimensional vectorial space. We show that this 
solution depends continuosly on q. Consider the linear operator 

-A A : (v 2 (n)), -»• (K, 2 (0)) A 

that associate to every u of (Vq(SI))\ the unique element — A\(u) such that: 

Vw G (Vo 2 (0)) A , / -A A ii v dx = [ -Au v dx. (11) 
Jn Jn 

So, — A A u is the orthogonal projection of —Am on (V r 2 (O))A- Observe that 
— Aa is a linear operator that acts on a finite dimensional vector space with 
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Ker(— A\) = {0}, so it is invertible. 

Now, let ei(x), e n (x) be an orthogonal base of (Vq(Q))\, and consider the 
function k : £1 — > (Vq(£1))\ that associates to every point 9 e O the unique 
function 6q t \ £ (Vo(£l))x defined as follows: 

n 

8 q ,\(x) = ^Jej(g)ej(x). 

i=l 

Observe that, by definition, Vv G (^o 2 (Q))a we have 

P P Tt Tl p 

/ 8 qt \vdx= I y^ei(q)ei{x)v{x) =y^ei(q) I ei(x)v(x) = v(q), 
Jn Jn l=1 i=1 Jo 

and, since u(q) = Jq 8 q (x)v(x) dx, we have 

Vv 6 (^ 2 (^))a, / VfrMaO dx= f S q (x)v*(x). (12) 

Let Uq t \{x) be a solution to[TU] Then, since — A>, is invertible and (|TTj) and (TT21 
hold, we have 

««,a(^) = A^ 1 o k(q). 

Since, as observed, k and (— A^) _1 are continuous functions, it follows that u q ^\ 
depends continuosly on q. Thus also 

F\(q) = E\(u g ,q) = - [ \\7u q ^{x)\ 2 dx + u{q) 



is continuous. 

Since f2 is compact, F\(q) has a minimizer which we denote by q\. 
Now let 



and 



q = lim q x 

AfU 



= limu g . a- 

AfU 



By Leibniz Principle, (u q ,q) is the minimizer of E(u,q) in £1 . Let us see that 
q £ £1* . By definition of Dirac ultrafunction we have that, for all q € <9f2, 5 q = 0, 
so u q (x) = and E(u q ,q) — 0, while E(u q ,q) < for every q 6 O*. So g € fi*. 
□ 



Remark 14 j4 similar problem that can be studied with the same technique is 
the problem of a electrically charged pointwise free particle in a box. Repre- 
senting the box with an open bounded set fl C R 3 , denoting by u q the electrical 
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potential generated by the particle placed in q E fl, then u q satisfies the Dirichlet 
problem 



The equilibrium point would be the point qo E £1 that minimizes the electrostatic 
energy which is given by 



namely, on the solution, the electrostatic energy is the opposite than the energy 
of a membrane-like problem in R 3 . In order to solve this problem we notice that, 
by definition of Dirac ultrafunction we have that, for all q E dfl, S q = 0. So 
Eei{q) > and E e i{q) = if and only if q E dtt. More precisely we have that 

• E e i{q) is infinite if the distance between q and dfl is larger than some 
positive real number; 

• E e i{q) is positive but not infinite for some q infinitely close to dfl; 

• E e i(q) = if and only if q £ dil. 
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